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Abstract: 

A measurement of the inclusive cross section for the deep-inelastic scattering of 
positrons off protons at HERA is presented at momentum transfers 8.5 < Q 2 < 
35 GeV 2 and large inelasticity y = 0.7, i.e. for the Bjorken-x range 0.00013 < 
x < 0.00055. Using a next-to-leading order QCD fit to the structure function F 2 
at lower y values, the contribution of F 2 to the measured cross section at high y is 
calculated and, by subtraction, the longitudinal structure function Ft, is determined 
for the first time with an average value of F L = 0.52 ± 0.03 (stat) 1 ;| 2 (syst) at 
Q 2 = 15.4 GeV 2 and x = 0.000243. 
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1 Introduction 



Precise measurements of the inclusive scattering cross section at the ep collider HERA are 
important for the understanding of proton substructure. In the one-photon exchange approx- 
imation, which is valid in the kinematic domain explored here, the deep inelastic scattering 
(DIS) cross section is given by the expression 

(J2{T 2lTnI [(2(1 -y)+y 2 )F 2 (x,Q 2 )-y 2 F L (x,Q 2 )}. (1) 



dxdQ 2 Q A x 

Here Q 2 is the squared four-momentum transfer, x denotes the Bjorken scaling variable, y = 
Q 2 / sx is the inelasticity, with s the ep center of mass energy squared, and a is the fine structure 
constant. The structure functions F 2 and Fl are related to the cross sections ot and <jl for 
the interaction of transversely and longitudinally polarized virtual photons with protons. In 
the Quark Parton Model F 2 is the sum of quark and antiquark distributions multiplied by x 
and weighted with the square of the electric charges of the quarks, while Fl is predicted to be 
zero for spin 1/2 partons [I]]. In Quantum Chromodynamics (QCD) Fl acquires a non zero 
value due to gluon radiation which is proportional to the strong coupling constant a s [0 with 
possibly sizeable higher order corrections in QCD perturbation theory ||. Measurements of 
Fl, expressed as the structure function ratio 

R=— ^— =°Jl (2) 
F 2 — Fl <Jt 

have been made by various fixed target lepton-hadron scattering experiments at higher x values 
[[]]. |5|. This paper presents the first determination of Fl(x,Q 2 ) at HERA in the deep inelastic 
region of 8.5 < Q 2 < 35 GeV 2 and very small x values between 1.3 ■ 10~ 4 and 5.5 ■ 10 -4 . 

The HI collaboration has recently reported a measurement of the structure function F 2 
in the range 3 ■ 10" 5 < x < 0.32 and 1.5 < Q 2 < 5000 GeV 2 , using data taken in the year 1994. 
The measurement was restricted to y values between 0.01 and 0.6 where the contribution of F 2 
to the cross section, Eq. ([!]), dominates. The F 2 values were extracted from the measured cross 
sections assuming theoretically computed values of Fl- A next-to- leading order (NLO) QCD 
analysis showed that the F 2 structure function can be well described by the DGLAP evolution 
equations @ in the kinematic range of the measurement. 

At high y the factors Y + = 2(1 — y) + y 2 and y 2 multiplying F 2 and Fl, respectively, are of 
comparable size. Therefore, the usual technique of extracting F 2 assuming a calculated F L is 
reversed and Fl is determined by subtraction of the F 2 contribution from the measured cross 
section. The following procedure is applied. Our measurement of F 2 ||, for y < 0.35, and fixed 
target data at larger x M are used to extract the parton distribution functions which are evolved 
in Q 2 according to the NLO DGLAP evolution equations. This provides predictions for the 
structure function F 2 in the high y region which allow, by subtraction of the contribution of F 2 
to the DIS cross section (cf. Eq. (1)), the determination of the longitudinal structure function 
Fl to be made. Note that the measurements of F 2 are well described by NLO QCD over four 
orders of magnitude in x and Q 2 while the evolution required here extends the maximum Q 2 
at fixed x by a factor of two only. Nevertheless, since an extended kinematic region is accessed 
here, where new effects could be important, it can not be excluded that the structure function 
F 2 behaves differently than assumed. 

Instead of subtracting the contribution of F 2 from the cross section one could perform a cross 
section analysis using the QCD predictions for both structure functions F 2 and Fl- This would 
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be conceptually different to the method employed in this paper because then an assumption 
would be required not only for F 2 but also for Ft which is less well known than F 2 . 

A salient feature of the subtraction method is a partial cancelation of systematic errors 
because the cross sections at low and at high y are measured using one common set of data. The 
experimental challenge is to measure the cross section at high y where the energy of the scattered 
positron E' e is comparatively low. The present measurement is made for 11 > E' e > 6.5 GeV, or 
0.6 < y < 0.78, which is an extension of the kinematic range covered by our previously published 
measurement of F 2 0. An understanding of the trigger efficiency, positron identification, 
photoproduction background and radiative corrections now becomes more demanding. 

The paper is organized as follows. Section 2 discusses the cross section measurement with 
particular emphasis on the high y region. Section 3 describes the QCD fit used to define the 
F 2 contribution for subtraction and presents the final results. A short summary is given in 
Section 4. 



2 Cross Section Measurement 
2.1 Kinematics 

In 1994 HERA was operated with positrons of energy E e = 27.5 GeV and protons of energy 
E p = 820 GeV. The event kinematics were reconstructed using the energy of the scattered 
positron E' e and the polar angle 9 e according to the relations 

Ql = f« = 1 " W (3) 

J- — He &e 

Here 8 e is defined with respect to the proton beam direction, defining the z axis, and x is 
calculated as Q 2 /sy with s = 4E e E p . At high y the determination of Q 2 and y from the 
reconstructed positron, rather than from the final state hadrons or a combination of both, is 
preferred because of the superior resolutions in Q 2 and x. The determination of the inclusive 
event kinematics using the variables E' e and 9 e is subsequently referred to as the "electron 
method". 

The previously published analysis || used the "sigma method" to determine F 2 (x, Q 2 ) for 
y < 0.15. This method combines the positron with the hadronic measurement by defining 

2 _ E' e 2 sm 2 9 e _ y h 

- — ; Vt, - — , (4J 

1 - Vt, l + y h -y e 



which avoids the resolution degradation of y e at low y. Here 

Ei(Ei-p z ,i) 

Vh = 2£ ' ( 5 ) 

where Ei and p Z} t are the energy and longitudinal momentum component of a particle i. The 
summation extends over all hadronic final state particles and the masses are neglected. 

The kinematic region of the Ft measurement was limited by the constraints 0.6 < y < 0.78, 
155° < 6 e < 171°. It was divided into six intervals of Q 2 with the limits (7.5, 10.0, 13.3, 
17.8, 23.7, 31.6, 42.2) GeV 2 and with central values chosen to be (8.5, 12.0, 15.0, 20.0, 25.0, 
35.0) GeV 2 or the corresponding values of x — Q 2 /sy at y — 0.7. 
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2.2 The HI Detector 



The HI detector || is a nearly hermetic apparatus built to investigate high-energy ep interac- 
tions at HERA. The measurement of the inclusive deep inelastic cross section relies essentially 
on the inner tracking chamber system and on the backward electromagnetic and liquid argon 
calorimeters. A superconducting solenoid surrounds both the tracking system and the liquid 
argon calorimeter, providing a uniform magnetic field of 1.15 T. 

The energy of the scattered positron was measured in the backward electromagnetic calori- 
meter (BEMC) behind which a scintillator hodoscope (TOF) was placed to veto proton beam 
induced background interactions. The identification of the scattered positron and the mea- 
surement of the polar angle made use of the backward multiwire proportional chamber (BPC) 
which was attached to the BEMC. In the kinematic range of this Fl measurement the positron 
angle was limited to 155° < 6 e < 171°. For these angles the scattered positron traversed the 
inner cylindrical proportional chamber (CIP) which could therefore be included in the positron 
identification requirement. The inner and outer proportional chambers (CIP at 18 cm radius 
and COP at 47 cm radius) were used in the trigger to reconstruct tracks of particles originating 
from the interaction region, and thus to reduce beam induced background events. The inter- 
action vertex was determined with the central drift chambers and the hadronic final state was 
reconstructed with the Liquid Argon calorimeter and the tracking detectors. 

The luminosity was determined from the cross section of the elastic bremsstrahlung process, 
ep — > ep7, measured with a precision of 1.5%. The integrated luminosity for this analysis is 
1.25 pb -1 . The final state positron and the photon scattered at very low Q 2 can be detected 
in calorimeters ("electron and photon taggers") which are situated 33 m and 103 m from the 
interaction point in the positron beam direction. 

The use of the HI detector for the inclusive DIS cross section measurement is further 
discussed in [||. A detailed technical description of the apparatus can be found in ||. 

2.3 Trigger 

The DIS event selection was based on events triggered in the BEMC by an energy deposition 
of more than 6 GeV, combined with a TOF requirement and a valid CIP-COP track signal. 
The efficiency of the BEMC energy requirement was monitored using a central track trigger 
and found to be better than 99% for the whole analysis region. The CIP-COP trigger required 
at least one track pointing to the interaction region. The efficiency was 96% after all selection 
cuts. It was found to vary little over the region of acceptance and to be well reproduced by 
the simulation of the trigger response. It was monitored down to a positron energy of 7.5 GeV 
by an independent BEMC trigger and was evaluated between 6.5 and 7.5 GeV by studying 
its dependence on 8 e , on the hadronic angle and on the charged track multiplicity comparing 
simulation with data. 

2.4 Event Selection 

The event selection criteria are summarized in Table |l|. The positron was identified as the most 
energetic cluster in the BEMC associated with a signal in the preceding BPC and, if geomet- 
rically accessible, in the CIP. For the determination of the event kinematics and background 
suppression a vertex had to be reconstructed with more than one track in the central drift 
chamber. 
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6.5 GeV< E' e < 11 GeV 
ei < 4 cm 
e 2 < 3.5 cm 
€3 < 5 cm 
—25 cm< z„tx < 35 cm 
N tr > 1 



Table 1: Summary of event selection criteria. For positron identification three estimators were 
used - 1\: reconstructed positron cluster radius in the BEMC; e 2 : distance from the closest BPC 
hit to the centroid of the positron cluster; e^: distance from the positron candidate trajectory 
to the closest active CIP pad (not used for tracks outside the CIP acceptance region). z vtx 
denotes the z position of the reconstructed interaction vertex and N tr is the number of charged 
tracks reconstructed in the central drift chambers. 



Deep inelastic events were generated using the DJANGO [l(J program which is based on 



HERACLES iTTfl for the electroweak interaction and on LEPTO [121 to simulate the hadronic 



final state. Photoproduction background was generated with the PHOJET [ 13 1 program. The 



detector response was simulated using a program based on GEANT |TJ| . The simulated events 
were subjected to the same reconstruction and analysis chain as the data. For comparisons with 
experimental distributions, all simulated spectra were normalized to the measured luminosity. 

Fig. la shows the distribution of the energy of the scattered positron for the events passing 
all selection criteria in the high y region. The experimental distribution is very well described 
by the superposition of the simulated spectra from DIS events and of the photoproduction 
background, discussed below. 

For the genuine DIS events at high y, the current jet particles are on average emitted 
backwards with respect to the proton beam direction. Thus there is a possibility that the 
largest energy cluster is not due to the scattered positron but to a hadronic energy deposition 
in the BEMC. In a Monte Carlo simulation of DIS events 3% of the selected positron candidates 
were found to be produced by the hadronic final state. In more than 99% of the simulated 
events the genuine positron was either the highest energy or second highest energy cluster. 
For a comparative study of the effect of misidentification at low energies, a positron finding 
algorithm was used which accepted an event even if the highest energy cluster failed to satisfy 
the selection conditions but the second highest energy cluster fulfilled them. The resulting cross 
section agreed to within 1% with that based on the standard positron finding algorithm which 
used the highest energy cluster. Fig. |b shows the BEMC energy distribution of the cluster 
with second highest energy which is well reproduced by the simulation. The background due to 
photoproduction is small because there is only a small probability to generate two high energy 
clusters in the BEMC in such events. 



2.5 Photoproduction Background 

At low energies and for large polar angles of the scattered positron there are two major sources 
of background in the candidate DIS events. Non-ep background occurs due to beam interac- 
tions with residual gas and beam-line elements. An effective filter against such events is the 
requirement of a reconstructed event vertex in the interaction region. The number of remaining 
beam-induced background events was estimated to be 2.5% in the lowest Q 2 interval and below 
1% everywhere else using non-colliding bunch events. 

The second, more difficult, background source is photoproduction, including low Q 2 DIS 
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events, in which the scattered positron escapes undetected along the beam pipe and in which 
an energy cluster from the final state particles fakes a positron signal in the BEMC. The 
typical characteristic of the 7p background is a rapid rise of the cross section towards lower 
cluster energy. Most frequently, the energy cluster in the BEMC is produced by a 7r° decay to 
two photons or by charged hadrons, mainly n ± . Energy clusters due to neutral particles are 
effectively removed by demanding a track pattern in the CIP and the requirement that a BPC 
hit coincide spatially with a BEMC cluster. Hadronic clusters typically have large cluster radii 
in the BEMC, and are rejected by the cut on e±, see Table 1. 

The Monte Carlo simulation was used to subtract bin by bin the remaining photoproduction 
background. A fraction of photoproduction interactions had the genuine final state positron 
detected ( "tagged" ) in the electron tagger. Fig. ||a shows the BEMC energy cluster distribution 
for such events which passed the selection criteria. Within the accepted range in energy (6.5 
to 11 GeV) the simulation reproduces the observed rate to within 4%. 

A further study was based on events with large energy cluster radius (ei) or without CIP 
validation (63). These samples predominantly consist of photoproduction events. The event 
sample rejected by the ei cut (Table 1) allows the study of faked positron signals by charged 
hadrons. The shape of the energy spectrum agrees well with the simulated distribution and the 
normalization is reproduced to within 7%. The event sample rejected by the €3 cut allows the 
study of 7T° induced background. The fake positron energy distribution of events in this sample 
is shown in Fig. §3. The normalization agrees to within 2% of the simulated rate. 

The photoproduction background amounts to < 20% for the lowest Q 2 interval and decreases 
to < 5% for the highest Q 2 interval. The normalization uncertainty was estimated to be 20% 
taking into account the fluctuations per bin of the simulated sample and also the fact that 
only part of the photoproduction events, with positron energies between 5 GeV and 15 GeV, 
is tagged. 



2.6 Cross Section Determination 

The deep inelastic scattering cross section was obtained by correcting the background sub- 
tracted number of events with the acceptance calculated from the Monte Carlo events, nor- 
malized to the measured luminosity. The cross section was corrected for higher order QED 
and electroweak contributions using the HECTOR []T3]| program. Starting from the GRV [|1(J 



parton distributions, a two step iterative analysis was performed to calculate the acceptance 
and the radiative corrections. This maintains the uncertainty of the cross section measurement 
due to input structure function variations below 1%. 

The radiative corrections were calculated to order a 2 with soft photon exponentiation 



TE, [17]]. Taking into account the hadronic track requirements they are about 35%. Detailed 



comparisons were made between the HECTOR result and the HERACLES [IT] Monte Carlo 
simulation which showed agreement at the per cent level. A study has also been made compar- 
ing the cross section results with and without a selection > 0.1 which, when applied, reduces 
the radiative corrections to about 15%. The resulting cross sections agreed to within 2%. 
The systematic error on the cross section is derived from the following contributions: 

• A 1% uncertainty of the BEMC energy scale leads to an error of about 1.5%. 

• A 1 mrad uncertainty of the measured polar angle of the positron causes a 2% error. 

• The radiative corrections lead to a cross section uncertainty of 2%. 
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• The vertex reconstruction efficiency is known to 2% apart from the lowest Q 2 interval 
where a 4% error was estimated. 

• As in U an efficiency error of 2% is assigned to account for global event selection, BPC 
efficiency and TOF veto uncertainties. 

• The various cut efficiencies, studied using different deep inelastic and background enriched 
data and simulated samples, lead to an estimated systematic error of 3% including the 
trigger efficiency error. 

• An extra error of 1% is estimated for positron misidentification effects using the Monte 
Carlo simulation and the study of the stability of the measurement against ignoring or 
considering the second highest energy cluster in the BEMC. 

• The photoproduction background was known to within 20%. This leads to a 4% error at 
Q 2 = 8.5 GeV 2 decreasing to 1% at Q 2 = 35 GeV 2 . 

Statistical errors in the Monte Carlo acceptance and efficiency calculations were computed and 
added quadratically to the systematic error. The total systematic error on the cross section is 
about 8% at Q 2 = 8.5 GeV 2 and about 6% for the higher Q 2 values. Most of the error sources 
scale with E' e and are only weakly dependent on 9 e in the range of the Fl determination. The 
statistical error is about three times smaller than the systematic error. 



Q 2 /GeV 2 


X 






^syst 


F 2 


A s tat 


^syst 


Rcalc 


8.5 


0.000135 


1.165 


0.027 


0.095 


1.354 


0.031 


0.110 


0.45 


12.0 


0.000190 


1.198 


0.026 


0.075 


1.375 


0.030 


0.086 


0.40 


15.0 


0.000238 


1.368 


0.032 


0.079 


1.561 


0.037 


0.090 


0.38 


20.0 


0.000317 


1.276 


0.034 


0.071 


1.445 


0.038 


0.080 


0.35 


25.0 


0.000396 


1.439 


0.042 


0.079 


1.651 


0.048 


0.091 


0.39 


35.0 


0.000554 


1.435 


0.062 


0.077 


1.634 


0.071 


0.088 


0.37 



Table 2: Inclusive cross section a = d 2 a/dxdQ 2 , eq.(l), scaled by the kinematic factor k = 
Q 4 x/ (27ra 2 • Y + ) with statistical and systematic errors. The Q 2 , x values correspond toy = 0.7 
in all bins. Also quoted are the values of F 2 corresponding to these cross section measurements 
with calculated R values, given in the rightmost column. The values of R = R ca ic were obtained 
using the GRV parton distributions [16] as input. There is an additional, overall normalization 
uncertainty of 1.5% due to the luminosity measurement error. 

The ep cross section is given in Table [| for the six new intervals in Q 2 and x. The analysis 
was also extended into the region of our previously published F 2 results for 0.6 > y > 0.03. 
In Fig. |3] the present cross section measurement is shown together with the data ||. There is 
everywhere good agreement in the region of overlap. The cross section is quoted in the form 
cr • Q 4 x/(27ra 2 • Y + ) = F 2 — y 2 F L /Y + which for small y is about equal to F 2 independently of 
Fl. The three lines drawn in Fig. 3 represent cross sections calculated using the QCD fit for 
F 2 which is described in the subsequent section and three different assumptions on Fl- The 
dashed-dotted and dashed lines correspond to the limits Fl = and Fl = F 2 , respectively, 
as required by the positivity of the cross sections ol and <tt- The solid line represents the 
cross section with F L calculated using the gluon and quark distributions obtained by the QCD 
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analysis of F 2 . It becomes apparent in Fig. 3 that at lowest x, corresponding to the high 
y region of this data, the cross section becomes very sensitive to the longitudinal structure 
function. On the contrary, at larger x, for about y < 0.35, the F 2 contribution dominates and 
the three lines nearly coincide. Most of the previously published F 2 data points are insensitive 
to the assumptions on Fl. 

In the publication || the measured DIS cross section was used to determine the structure 
function F 2 assuming R to be given by the GRV parton distributions [15 using the relation JTPI . 
This measurement represents an extension of the previous cross section data towards lower x. 
In order to provide a consistent set of structure function values the corresponding six values of 
F 2 were derived following the same procedure (see Table 2). Note that these values are rather 
sensitive to the R values chosen. 



3 Determination of Fl 
3.1 QCD Fit 

For the subtraction of the F 2 contribution to the cross section a NLO QCD fit was performed 
using the DGLAP evolution equations. The fit used the HI data || for y < 0.35. The BCDMS 
proton and deuterium data || were used to constrain the high x behaviour of the parton 
distributions. In contrast with the previous QCD analysis performed by HI B, the NMC data 
pO were not included in the standard fit to ensure a maximum weight of the HI data in the fit 



procedure. The starting point of the evolution was chosen to be Ql = 5 GeV 2 and all data with 
Q 2 > Qmin = 1-5 GeV 2 were included in the fit. To avoid possible higher twist effects, BCDMS 
data in the range x > 0.5 for Q 2 < 15 GeV 2 were not included in the fit. The normalization 
of the HI data was kept fixed. The fit used three light flavors with the charm contribution 



added using the NLO calculation of the photon-gluon fusion process [jig , |21fl . Furthermore, the 
momentum sum rule was imposed and the integral over the valence quark distributions was 
set to 2 for u v and to 1 for d v . The input parton distributions at the starting scale Qq were 
parameterized as follows: 



xg(x) 


= A g X B9 ( 


1 


— x] 


° 9 , 






= A u x Bu 


(1 


— X 




[l + D u x + E u Vx), 


ccdu (^cc j 


= A d x Bd ( 


1 


— x) 




1 + D d x + E dv ^), 


xS(x) 


= A s x Bs 


(1 


— X 




(1 + Dsx + EsV^), 



(6) 

where S = u = d = 2s defines the sea distributions. Three fits with different, fixed Aqcd values 
were performed. The best -y^/ndf of 506/(505-15) was obtained for Aqcd = 210 MeV. The 
fitted parton distribution functions were evolved into the new domain using the NLO DGLAP 
equations and used to calculate the corresponding values of F 2 . 

Table [3] summarizes the Q 2 averaged uncertainties in F 2 arising from the fit procedure. The 
total uncertainty due to the fit assumptions amounts to 1.7%. The resulting absolute error of 
the longitudinal structure function AFl ~Y + /y 2 - AF 2 is approximately 0.07. 

There is a small dependence of the structure function F 2 on Fl due to the assumption made 
for Fl in the cross section analysis for y < 0.35. Thus the two extreme assumptions Fl — 
and Fl = F 2 were used and two modified structure functions F 2 were derived as input to two 
QCD fits. This changed the QCD predicted F 2 at y — 0.7 on average by —1.6% and +3.8%, 
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fit assumption 


uncertainty in % 


NMC data used 


1.4 


change of Aq C d by 50 MeV 


0.7 


g(x,Ql)-(l + E^) 
Q 2 min = 5 GeV 2 
Ql = 3 GeV 2 


0.1 
0.6 
0.4 



Table 3: Uncertainty, relative to the result of the standard fit (see text), of the structure 
function F 2 averaged over the Q 2 range of the Fl data for various assumptions in the QCD fit 
procedure. 



respectively. Thus an asymmetric error on F L was introduced and was added in quadrature to 
the other systematic errors. 

Two different cross checks were made of the prediction of F2 at lowest x by using the 



perturbative dipole model with kr factorization [22] and an empirical model based on the 



similarity of the rise of F2 at low x and the evolution of the charged multiplicity with energy 



in e + e~ collisions p3|| . The model parameters were determined using the previously published 
HI F 2 data || for y < 0.35 and Q 2 boundaries given by the limitations of these approaches. 
The three-parameter F 2 function in the dipole model, calculated at y = 0.7, is only 2% lower 
than the structure function obtained by the evolution procedure described above. Similarly, 
the two-parameter F 2 function of the empirical model is on average 2% higher than the QCD 
fit result at y = 0.7. Thus both approaches to extrapolate F 2 would lead to a result for the 
longitudinal structure function in very good agreement with the one obtained subsequently 
using the QCD fit for the description of F 2 . 



3.2 Results 

The measured longitudinal structure function Fl is given in Table f|. The systematic error 
consists of the following contributions: 

• The experimental errors of the cross section measurement which are uncorrelated with 
the error of the data entering the QCD fit at lower y. These are error sources, such as 
the tracking trigger or CIP efficiency, which are specific for the high y range and to this 
analysis. 

• The error due to possible variations of the assumptions in the QCD fit procedure as 
discussed in Section 3.1. 

• The experimental errors like energy and angle uncertainties and global efficiency and 
luminosity errors which are mostly common to both the low and the high y region. A 
correlation of these errors is introduced through the QCD fit of F 2 . This error source 
includes also the statistical error of the fit result. 

The third contribution comprises several effects. For example, any global shift common to 
the high y data and the HI data used in the fit, like the luminosity uncertainty, gets reduced 
to about 1/3 of its magnitude. A reduction of the error is observed as well for the polar angle 
uncertainty. However, the error of the energy of the scattered positron is not compensated. 
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The HI F 2 data || for y > 0.15 were obtained with the electron method and those at smaller 
y with the sigma method. A 1% increase in E' e increases the cross section for y > 0.15 and 
decreases it below that value, a behaviour which leads to a large change in the \ 2 °f the fit. 
Thus the measurement of the positron energy is the dominating "correlated" error although it 
has only a small effect on the cross section at high y. 



Q 2 /GeV 2 


X 


F L 


A stat 


~~^~A S y S t 


A S y S t 


AUTIC 

exp 


A cor 
exp 




8.5 


0.000135 


0.51 


0.06 


0.29 


0.27 


0.17 


0.19 


0.06 


12.0 


0.000190 


0.63 


0.06 


0.28 


0.25 


0.15 


0.18 


0.07 


15.0 


0.000238 


0.35 


0.08 


0.29 


0.27 


0.15 


0.19 


0.08 


20.0 


0.000317 


0.67 


0.08 


0.28 


0.26 


0.14 


0.18 


0.07 


25.0 


0.000396 


0.33 


0.10 


0.25 


0.22 


0.15 


0.14 


0.07 


35.0 


0.000554 


0.39 


0.15 


0.24 


0.21 


0.14 


0.14 


0.07 



Table 4: The longitudinal structure function Fl with statistical (A sta t ) and systematic errors 
(±A syst ): A"™p is the uncorrected experimental cross section error at high y, A c °^ p is the 
correlated experimental error and A fu is the error introduced by the QCD fit uncertainty The 
total systematic error contains also the asymmetric contribution due to the assumptions on Fl 
in the determination of F 2 , see Section 3.1. The Q 2 , x values correspond to y = 0.7 in all bins. 

The six measurement values enable a determination to be made of the mean Fl and its 
derivative dF^/d In (x) for Q 2 = 0.7-sx from a straight line Fl = a+b-\n(x). Taking into account 
the error correlations between the six data points we obtain a mean Fl = 0.52 ±0.03 (stat) 1q]22 
(syst) and a derivative dF L /d\n(x) = -0.085 ± 0.080 (stat) tg;gjg (syst) at Q 2 = 15.4 GeV 2 
and x = 0.000243. Note that the derivative has comparable statistical and systematic errors 
while the error of the mean Fl is dominated by systematics. Fig. f| shows the data of Table 
[| and the extreme limits of Fl = and Fl = F 2 using the QCD fit. Without utilizing the 
measured dependence on x or Q 2 , these extremes are excluded with 2.3 and 4.0 times the total 
error, respectively. 

At low x the longitudinal structure function is related to the gluon distribution. The 
dashed band in Fig. 4 represents the calculation of Fl according to @ for three light quarks 



and according to ]21j for the charm contribution. The input gluon and quark distributions are 
determined by the NLO QCD fit described in Section 3.1. The width of this band is determined 
by the experimental errors of the F 2 data, taking into account their point-to-point correlations, 
and by the fit uncertainties discussed above. At the present level of accuracy there is consistency 
between the structure function Fl determined from this analysis and that calculated from the 
gluon and quark distributions. 



4 Summary 

Based on data taken in 1994 with a luminosity of 1.25 pb -1 , an inclusive measurement of the 
deep inelastic cross section measurement at y = 0.7 has been used to determine for the first time 
the longitudinal structure function Fl(x, Q 2 ) at very low Bjorken x. The analysis assumed the 
proton structure function F 2 (x, Q 2 ) to be in accordance with next-to-leading order perturbative 
QCD. The result excludes the extreme limits of Fl = and Fl = F 2 , corresponding to R = 
and R = oo, by 2.3 and 4.0 times the total error on Fl- The result is consistent with a higher 
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order QCD calculation of which essentially relied on the gluon distribution as determined 
from the F 2 structure function data. 
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Figure 1: Energy distributions of a) the highest energy and b) the next highest energy BEMC 
clusters for the final data sample. The simulated spectra are normalized to the luminosity of 
the data. 




Figure 2: Energy distributions of the highest energy cluster in the BEMC a) for photoproduc- 
tion events in which the scattered positron was tagged and b) for the events rejected by the 
CIP requirement. The simulated spectra are normalized to the luminosity of the data. 
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Figure 3: Double differential cross section ndu / dxdQ 2 = F 2 — y^Fl with k = Q 4 x/(2ira 2 ■ Y + ) 
in six Q 2 bins as a function of x. For y > 0.6 this analysis (open points) extends the previously 
published measurement [6] (closed points) towards lower x and is drawn here with full errors. 
The open points at larger x are given without errors for ease of comparison with the data of [6]. 
The three lines represent calculated cross sections using for F 2 the QCD fit, as described in sect. 
3.1, and three different assumptions for F^. These are the two extremes, F^ = (dashed-dotted 
line) and Fi = F 2 using F 2 from the QCD fit (dashed line), and Fi as calculated in NLO from 
the quark and gluon distributions determined by the QCD fit (solid line). 



16 



2 

F L 

1 .5 









10 




20 


Q 2 


/GeV 2 


I II 

: H1 








1 


1 r 




< 




II 


> 






< 


























» 














i i i i i i 



10 4 3-10 



Figure 4: Longitudinal structure function Fl as function of Q 2 or x = Q 2 /sy for y = 0.7. 
The inner error bars are the statistical errors. The full error bars represent the statistical 
and systematic errors added in quadrature. The error band represents the uncertainty of the 
calculation of Fl using the gluon and quark distributions, as determined from the NLO QCD 
analysis of the HI data [6] for y < 0.35 and the BCDMS data [8]. The dashed lines define the 
allowed range of Fl values from Fl = to Fl = Fi where F2 is given by the QCD fit. 
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